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We analyze optical spectroscopy data of the electron-doped superconductor (Pr2- a; Ce a; )Cu04 
(PCCO) to investigate the coupling of the charge carriers to bosonic modes. The method of anal- 
ysis is the inversion of the optical scattering rate T~p(uj,T) at different temperatures T by means 
of maximum entropy technique combined with Eliashberg theory. We find that in the supercon- 
ducting state the charge carriers couple to two dominant modes one at ~ f 2 meV and a second 
one at ~ 45 meV as well as to a high energy background. The low energy mode shows a strong 
temperature dependence and disappears at or slightly above the critical temperature T c . The high 
energy mode exists above T c and moves towards higher energies with increasing temperatures. It be- 
comes less prominent at temperatures > 100 K above which it evolves into a typical spin-fluctuation 
background. In contrast to the hole-doped High-T c superconductors PCCO proves to be a super- 
conductor close to the dirty limit. 

PACS numbers: 74.20.Mn 74.25.Gz 74.72.-h 



I. INTRODUCTION 



The optical scattering rate t~ v (w) in the infrared 
regime evolved into an important tool to extract micro- 
scopic information on the coupling of charge carriers to 
bosonic modes in the form of the electron-exchange bo- 
son interaction spectral density I 2 x{u) in the High-T c 
cuprates. The definition of the optical scattering rate 
itself is based on a generalized Drude form valid for cor- 
related electron systems: 
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This equation relates the complex optical self energy 
E op (o;,r) at a given temperature T to the complex op- 
tical conductivity a{w,T). The optical scattering rate is 
related to the imaginary part of the optical self energy 
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Here Q p is the plasma frequency and u> an energy vari- 
able. 

It was first demonstrated by Marsiglio et al. 1 that 
there exists an approximate relation between I 2 x(oj) and 
T~p (uj) of the form: 



(3) 



This relation is valid only at low temperatures and up 
to energies at which W(u>) first becomes negative, i.e. 



unphysical. This second derivative method was used by 
Carbotte et al. 2 to demonstrate that in the optimally 
doped system YBa2Cu306+5 (YBCO) the quasiparticles 
couple to a boson resonance at 41 meV which corresponds 
to a spin one resonance observed by neutron scattering 
in the imaginary part of the spin susceptibility. 3 ' 4 More- 
over, the temperature dependence of this resonance peak 
was found to be identical to the T variation of the equiv- 
alent structure in the spectral function J 2 x(w) derived 
from optics. 5 Finally, Schachinger and Carbotte 6 pre- 
dicted such a spin one resonance to exist in the thallium 
compound T^P^CuOg+a (T12201). This was later con- 
firmed by He et al. 7 using neutron scattering. 

The next step in the development of methods which 
help to extract information on the electron-exchange bo- 
son interaction spectral density I 2 x(oj) from optics, was 
made by Dordevic et al. s who developed a new method 
based on inverse theory. These authors concentrated on 
the approximate relation 
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reported by Shulga et al. 9 Eq. (4) is based on Eliashberg 
theory and is valid in the normal state. Here Tj" is 
an energy independent impurity scattering rate and the 
kernel K (uj, v\ T) is given by: 
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Relation (4) was extended by Carbotte and 
Schachingcr 10 to the superconducting state of a 
d-wave superconductor at T — using the kernel 

K{u,v;T = Q) = — ((w-i/)0[w-2A(0) - v] 




This kernel is based on a clean limit, i.e.: r imp = 0, per- 
turbation theory expansion of BCS theory reported by 
P.B. Allen 11 for an s-wave superconductor which has also 
been considered by Dordevic et al. 8 In as much as one can 
think of a d-wave superconductor as a superposition of 
s-wave with variable gaps, Eq. (6) follows as a first ap- 
proximate generalization of Allen's work to <i-wave. Here 
{■■■}■& denotes the average which can be limited to the 
interval •& £ [0, 7r/4] for symmetry reasons. Furthermore, 
A($) = A cos(2$) reflecting the d-w&ve symmetry of the 
superconducting order parameter. Eq. (6) ensures that 
the optical scattering rate is finite in the superconduct- 
ing state for uo > 0. Finally, E(x) is the complete elliptic 
integral of second kind and 9(x) is the step function. 

Eq. (6) explains why in the superconducting state of a 
d-wave superconductor the signatures of the I 2 x(y) spec- 
tral density in the optical scattering rate t~^(oj) are dis- 
placed at least by the gap amplitude A. 12-14 This is in 
contrast to a classical s-wave superconductor in which 
those signatures are displaced by 2A , with A the gap- 
edge. Finally, we would like to point out that according 
to Eq. (4) non-zero contributions to the positive definite 
spectral density I 2 x{v), i.e.: non-zero coupling of charge 
carriers to an exchange boson at energy v, will always re- 
sult in an increase of the optical scattering rate because 
both kernels, Eqs. (5) and (6), are positive definite. Thus, 
depressions or pronounced peaks in the optical scatter- 
ing rate r" 1 (to) cannot be explained by electron-exchange 
boson interaction and indicate additional influences like, 
for instance, infrared activated phonons or quasiparticle 
density of states effects. 

The spectra derived from the deconvolution of Eq. (4) 
still contained negative and, thus, unphysical parts. Nev- 
ertheless, it was demonstrated by Dordevic et al. 8 that 
the application of a deterministic constraint can result 
in the regularization of the problem which allows to the 
removal of unphysical negative values from the inverted 
spectral function l 2 x{v). The authors also demonstrated 
that such a regularization can result in a reduced quality 
of data reproduction. 

In a final step, so far, in the development of inversion 
techniques it was demonstrated by Schachinger et al. 15 
that the application of maximum entropy techniques 16 
maximizing the generalized Shannon- Jaynes entropy 17 
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to deconvolute Eq. (4) will result in a positive definite 
spectral density I 2 x(v) ensuring best possible data re- 
production. Here, the entropy S measures the distance 
of a candidate vector a = {aj\j = 1, . . . , N} from the de- 
fault model vector m = {rrij\j = 1, . . . , N} which repre- 
sents the most probable solution prior to the observation 
of any data. In case of insufficient background informa- 
tion it should be chosen constant, i.e.: mj = const, Vj, 
as will be done throughout this paper. 

To compensate for the approximate nature of Eq. (4) 
an additional least squares fit procedure based on Eliash- 
berg equations which have been extended to <i-wave su- 
perconductors (see Appendix) proved, finally, to be very 
successful in the inversion of optical Bi 2 Sr 2 CaCu208+(5 
(Bi2212) data at various temperatures and doping 
levels. 18 

All this research resulted in one common denominator: 
in YBCO, Bi2212, and T12201 the optical data suggest 
coupling of the charge carriers to a pronounced boson 
mode, an 'optical' resonance, at energies which agree 
in most cases with the energies at which a spin one 
resonance is found by neutron scattering in the imagi- 
nary part of the spin susceptibility. In a new experi- 
ment Vignolle et al. 19 showed, using neutron scattering, 
that the imaginary part of the spin susceptibility in opti- 
mally doped (La2-xSr x )Cu04 (LSCO) develops two peak 
structures, a resonance at the low energy of 12 meV while 
the second peak was found at <~ 50meV. A maximum 
entropy analysis of optical data reported by Gao et al. 20 
on epitaxial optimally doped LSCO thin films confirmed 
that, indeed, the electron-boson spectral density I 2 x{u) 
develops two peaks at the energies reported by neutron 
scattering. 21 This is also in agreement with earlier results 
suggested by Zhou et al. 22 from angular resolved photo 
emission spectroscopy (ARPES) experiments in under- 
doped samples. They interpreted their results as a cou- 
pling of the charge carriers to phonons which is in some 
contrast to the results reported by Vignolle et al. 19 and 
Gao et al. 20 

Neutron scattering experiments by Wilson et al. 23 
report the existence of a spin one resonance in 
the imaginary part of the spin susceptibility of 
(Pri_ a: LaCe a: )Cu04-5 (PLCCO) for x = 0.12 with its 
peak at ~ 11 meV. It is centered around (ir/2, ir/2) in the 
two dimensional CuO Brillouin zone as is the case in all 
other High-T c superconductors and it disappears at tem- 
peratures above T c . The doping (x) dependence of this 
peak was studied by M. Fujita et al. 24 Scanning tunneling 
microscope (STM) experiments performed on the same 
material by Niestemski et al. 25 support the existence of 
such a peak at <~ 10 meV. As optical data are available in 
the similar system (Pr 2 _ a; Ce a ;)Cu04_5 (PCCO) we con- 
centrate in this paper on the task to extract information 
on the electron-boson spectral density I 2 x(^) from op- 
tical scattering rates measured by Homes et al. 26 on a 
PCCO single crystal with x = 0.15 (T c = 20 K) and from 
similar data reported by Zimmers et al. 27 ' 28 on thin, epi- 
taxially grown PCCO films with x = 0.15 (T c = 21 K). 
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In Sec. II we discuss the samples used for the analysis 
in their main features. Data extraction techniques will 
also be discussed shortly within this section. Section III 
concentrates on the analysis of the optical data and the 
discussion of the important features of the I 2 x(w) spectra 
derived from optical data. Finally, Sec. IV gives a short 
summary. An appendix has been included to provide 
the necessary background information on the extended 
d-wave Eliashberg formalism used for impure systems. 



II. MATERIALS 



The optimally doped (x 
crystal 



0.15) PCCO single 



We concentrate first on the experimental results re- 
ported by Homes et al. 26 for an optimally doped PCCO 
single crystal (x = 0.15, T c ~ 20 K). The authors report 
two plasma frequencies, the first value tt p = 13 000 ± 
200 cm -1 (1.64 eV) reflects only those carriers which par- 
ticipate in coherent "Drude" transport, and a larger value 
for fi p = 19 300cm -1 (2.4 eV), which was determined 
from a modification of the finite-energy /-sum rule. 29 ' 30 
The two values of O p obtained using these different ap- 
proaches indicate that f2 p is not well known in this ma- 
terial, and that its value depends heavily on the method 
chosen to extract it from the experiment. This leads 
to some uncertainty because according to Eq. (2) the 
plasma frequency sets the scale of the optical scatter- 
ing rate. The smaller value of f2 p finds its justification 
in a two-component system where the Drude response is 
not directly connected to the mid-infrared region which is 
conceived to be due to another band of electrons; this was 
the approach used in the original analysis of the PCCO 
single crystal data. 26 However, in this work we use a one- 
component approach and as a result the larger value of £l p 
is more appropriate. In this view Drude and infrared re- 
gions come from the same electrons with the infrared part 
coming from the incoherent boson-assisted processes, and 
the Drude originates from the coherent quasiparticle re- 
sponse part of the carrier spectral function. Following 
this argument, the experimental r~p(uj) data has been 
derived from the raw reflectance data using il p = 2.4 eV 
together with the dielectric constant at infinity, eoo, set 
equal to four. It is worth remarking at this point that the 
dielectric function, e(w), and the conductivity are related 
through the expression (t(uj) — — iuj[e(uu) — e^/Air. Thus, 
too is required to derive the imaginary part of the optical 
conductivity, and as a consequence the optical scattering 
rate is also influenced by this important parameter. 

The optical scattering rate at the lowest normal state 
temperature reported (30 K) has a zero frequency off- 
set of ~ 22meV (~ 177 cm -1 ) which indicates substan- 
tial contributions from impurity scattering. In addition, 
there are sharp features in the scattering rate shown 
in Fig. 1(a) at approximately 306, 338, and 434 cm -1 
(« 37.9, 41.9, and 53.8 meV) which have been identi- 
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FIG. 1: (Color online) The optical scattering rate t„ v (u) in 
meV in an optimally-doped (x = 0.15) thin PCCO single 
crystal (T c = 20 K) for several temperatures above and below 
T c with (a) the three infrared-active E u modes at ~ 306, 338, 
and 434 cm" 1 (f» 37.9, 41.9, and 53.8 meV) present, and (b) 
with the same phonons removed. 



fied as the infrared-active E u lattice vibrations in this 
material. 26 > 31 ~ 33 In order to remove this sharp structure 
from the scattering rate, the dielectric function has been 
modeled using a series of Lorentz oscillators 



e(w) c 



E 



P,3 



(8) 



where ujj , jj and u> p j are the frequency, width and effec- 
tive plasma frequency of the j'th vibration, respectively. 
The oscillator parameters have been determined from fits 
to the real part of the optical conductivity using a sim- 
ple linear background; the fitted oscillator values are then 
used to generate the phonon contribution to the complex 
dielectric function which is in turn subtracted from the 
experimentally determined values for the real and imag- 
inary parts of the dielectric function, allowing the sharp 
phonon features to be removed from the optical scatter- 
ing rate, as shown in Fig. 1(b). 
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B. The optimally doped (x = 0.15) PCCO thin film 

We also analyzed data published on thin films by Zim- 
mers et al. 27,28 and concentrate on the x = 0.15 sam- 
ple. The reflectivity of a thin film also contains infor- 
mation on the optical properties of the underlying sub- 
strate. Hence one cannot obtain the film optical func- 
tions using a straight-forward Kramers-Kronig inversion. 
In the case of the PCCO films, the substrate response 
was determined separately. The reflectivity of the total 
(film plus substrate) system was then fitted by construct- 
ing a suitable Drude-Lorentz dielectric function for the 
film. 34 One should note that the parameters used in this 
Drude-Lorentz dielectric function are not to be taken as 
representative of intrinsic excitations. They should be 
looked at as a convenient parameterization of the ma- 
terial dielectric function. Once this dielectric function 
is known, any other optical function can be generated 
straightforwardly. 

Figure 2(a) presents the optical scattering rate t~^(uj) 
vs to for the optimally doped, epitaxially grown thin 
PCCO film (x = 0.15, T c = 21 K). The plasma frequency, 
calculated from a finite sum rule, is fl p — 17 570 cm -1 
(2.18 eV) and eoo was set to four. Figure 2(b) shows the 
scattering rate once phonons have been eliminated. 

The zero frequency offset of the optical scattering rate 
^(0) = 32.5meV'(- 262 cm- 1 ) at T = 25 K which 
suggests a much higher impurity concentration in com- 
parison to the single crystal discussed in the previous 
subsection. The overall shape of the scattering rate in 
the x = 0.15 film is similar to the crystal response. The 
phonons in the film tend to be broader and less screened 
by free carriers than in the crystal. These effects are con- 
sistent with a larger disorder and smaller weight of the 
coherent peak in the film, as we describe below. 

There is a scaling by roughly a factor of two required 
if one wants to compare directly the x = 0.15 film and 
crystal data. It comes from the difference in the ab- 
solute value of their respective scattering rates. It is 
worth noticing that this is not an effect originating in 
the scattering rate of carriers participating in the coher- 
ent transport. As discussed in Sec. II A, taking into ac- 
count the coherent transport alone, the low frequency 
optical conductivity of the x — 0.15 single crystal at 
T = 30 K can be described by a Drude peak having 
flp = 1.61 eV (13 000cm- 1 ) and t^~(u = 0) rj llmeV 
(90 cm- 1 ). Performing the same analysis on the x = 0.15 
film at T = 25 K one obtains fl p = 1.18 eV (9 500 cm" 1 ) 
and T op 1 (tj = 0) = llmeV. The scattering rate is the 
same but the weight of the Drude peak in the film is 
significantly smaller. As the film and the crystal com- 
positions are nominally the same, the /-sum rule states 
that the area under their respective optical conductivi- 
ties should be the same. Hence, the weight lost in the 
coherent peak of the film is redistributed as an incoher- 
ent background. This effect will produce, as observed, a 
higher mid-infrared optical conductivity in the film and 
hence a broader incoherent scattering rate. Thus, a fac- 
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FIG. 2: (Color online) The optical scattering rate t~^{w) 
im meV in an optimally doped (x = 0.15) PCCO thin film 
(T c — 21 K) for several temperatures above and below T c 
with (a) the infrared active phonon modes present and (b) 
with these phonons removed. 

tor of two is not far from the difference in dc- conductivity 
\a\(uj) extrapolated to u) = 0] observed in these two sam- 
ples (3xl0 4 f7- 1 cm- 1 for the crystal; 1.7 X lO 4 ^- 1 cm" 1 
for the film). 

III. DATA ANALYSIS AND DISCUSSION 
A. The x = 0.15 PCCO single crystal 

Because of the rather large normal state zero frequency 
offset of the optical scattering rate PCCO can no longer 
be treated in the clean limit as has been done so far for 
the systems Bi2212, LSCO, T12201, and YBCO. Thus, 
the impurity scattering rate r^p gains importance in 
Eq. (4) and has to be treated as an external parameter 
in the maximum entropy deconvolution of this equation. 

The complete inversion procedure consists of two steps. 
First Eq. (4) together with the appropriate kernel is de- 
convoluted using a classical maximum entropy method. 
The default model is set to a constant which is adjusted 
from temperature to temperature to ensure best possible 
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data reproduction at high energies u > 250 meV. This 
results in a first approximation for the electron-exchange 
boson interaction spectral density I 2 x(w)- In a second 
step this approximate spectral density is refined using a 
least squares fit procedure based on the full non-linear 
Eliashberg equations (A.l). This second step is of par- 
ticular importance when inverting superconducting state 
data. 

Fig. 3 reports the results of our data analysis. 
Figure 3(a) demonstrates the quality of data reproduc- 
tion when the I 2 x{^) spectra shown in Fig. 3(b) are 
used to calculate the optical scattering rate using the 
Eliashberg equations (A.l). The heavy lines represent 
the data and the thin lines theory. For the normal state, 
T = 295 K, 200 K, 100 K, and 30 K Eliashberg equation 
(A. lb) with K{v + i0 + ;z9) = was used to calculate 
the renormalized quasiparticle energies uj(v + i0 + ). In 
this case the impurity scattering rate r^ p is only a term 
which is added to the renormalized energies. An im- 
purity parameter t + — 3.1 meV (see Appendix) was re- 
quired for best over all data reproduction. It corresponds 
to an impurity scattering rate = 2nt + w 19.5 meV 
for all temperatures. This is in excellent agreement with 
the zero frequency offset of ~ 22meV reported for the 
T = 30 K data. 

Particular attention is required for the superconduct- 
ing state. First of all, impurities are always pair-breaking 
in d-wave superconductors and thus effectively reduce the 
critical temperature of the impure sample in comparison 
to the 'clean limit' critical temperature. Furthermore, 
there are two limits of impurity scattering to be con- 
sidered. One limit is described by unitary or resonant 
scattering which is characterized by the parameter T + in 
Eq. (A. lb). The other limit is Born's scattering (or weak 
scattering) characterized by the impurity parameter t + . 
In reality the scattering law is intermediate between uni- 
tary and Born limit scattering. This is characterized by 
the parameter c in Eq. (A. lb); c = is the unitary limit 
and c — > oo the Born limit. Another complication arises 
from the fact that the relevant kernel (6) does not con- 
tain impurity scattering in contrast to its normal state 
counterpart Eq. (5). Furthermore, Eq. (6) requires some 
knowledge of the size of the gap amplitude A. Fortu- 
nately, this is not critical in this particular case because 
the rather big impurity scattering rate of ~ 22 meV, sug- 
gests that the sample will already be in the gapless regime 
even at T — 5 K. 35 Thus, it will be sufficient for the first 
step, the maximum entropy deconvolution of Eq. (4), to 
use that value of A for which the best data reproduc- 
tion can be achieved. This procedure results in a first 
approximation of the I 2 x(u) spectrum which is then pa- 
rameterized and improved by a least squares fit to the 
data using the solutions of the full non-linear Eliashberg 
equations to calculate the optical scattering rate. In this 
step the impurity parameter t + stays unchanged and only 
Born scattering is treated but is included in each itera- 
tion of the Eliashberg equations (A.l) and, of course, this 
greatly reduces the value of the critical temperature T c . 
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FIG. 3: (Color online) (a): The optical scattering rate t^(uj) 
in meV as a function of u> (in meV) at five temperatures, 
namely T = 300 K (magenta, dash-double dotted line) , 200 K 
(blue, dash-dotted line), 100 K (green, dotted line), and 30 K 
(red, dashed line) in the normal state and T = 5K [black 
solid line, the only curves with t^(uj = 0) = 0] in the 
superconducting state for the optimally doped PCCO sin- 
gle crystal (x = 0.15). The heavy lines present experimen- 
tal data while the light lines correspond to theoretical re- 
sults obtained from solutions of the full Eliashberg equations 
(A.l) using the I 2 x(&) spectral densities shown in the bot- 
tom frame of this figure. The impurity parameter was set to 
t + — r + = 3.1 meV. (b): The electron-exchange boson inter- 
action spectral density I 2 x(^>) obtained from the inversion of 
the experimental data presented by the heavy lines in the top 
frame of this figure. 



The result are the spectral density I 2 x{<jj) shown by the 
solid line in Fig. 3(b). 

In a final step the low energy regime < cu < 10 meV 
in which the initial slope of (cj) is dominated by impu- 
rity scattering is to be fitted using full non-linear Eliash- 
berg theory, Eqs. (A.l), taking into account the different 
laws of impurity scattering. As the scattering rate itself 
has already been determined, the appropriate scattering 
law can now be determined by the best possible fit to 
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FIG. 4: (Color online) The solid black line represents the 
experimental t~^(u)) as a function of u) at T = 5K. The 
dashed curve presents the results of an Eliashberg theory cal- 
culation using unitary impurity scattering described by the 
impurity parameter r + = 3.15 meV. The dotted line cor- 
responds to Born limit scattering described by the param- 
eter t + = 3.15 meV. Finally, the dash-dotted curve repre- 
sents an intermediate case with T + = 3.94 meV and c = 0.5 
and, furthermore, the dash-double dotted curve belongs to 
r+ = 6.3 meV and c= 1. 



the data in this low energy region. Fig. 4 demonstrates 
the results of such a procedure. The solid black line 
corresponds to the experimental data, the dashed line 
represents unitary limit scattering with T + = 3.1 meV 
and the dotted line shows the result for Born limit scat- 
tering with t + = 3.1 meV. For comparison we include 
additional results for intermediate impurity scattering, 
namely T+ = 3.875 meV with c = 0.5 and T+ = 6.2 meV 
and c = 1. Obviously, impurity scattering is much closer 
to unitary than to Born type scattering but it is impos- 
sible to decide whether a better fit is found for c = 
(unitary scattering) or c = 0.5 (intermediate scattering). 
It also becomes apparent that impurity scattering does 
not affect the energy dependence of (uj) for energies 
uj > 30 meV and, hence, there is no need to get yet a new 
estimate for I 2 x(uj)- It is certainly interesting to note in 
passing that the clean limit (intrinsic) critical tempera- 
ture T c o for t + — is approximately 61 K so that the 
presence of impurity scattering has greatly reduced the 
value of the critical temperature. 

We return to Fig. 3(b) which presents the spectral den- 
sities I 2 x{uj) as a result of the inversion process. Five 
temperatures are shown, namely T = 5K (solid line), 
30 K (dashed line), 100 K (dotted line), 200 K (dash- 
dotted lines), and 295 K (dash-double dotted line). The 
main feature is a very pronounced peak centered around 
~ 45 meV and which gets reduced in amplitude with in- 
creasing temperature. This main peak is followed by 
a valley-hump structure which also becomes less pro- 



FIG. 5: (Color online) The single crystal experimental op- 
tical scattering rate Top (uj) as a function of uj at T = 30 K 
(black solid line). The two straight lines (red, dashed and 
solid) are used to emphasize the change in slope in r^, (ui) 
at uj — 45 meV. The (blue) dash-dotted line corresponds to 
the spectral density I 2 x(uj) derived from the r ( ^ ) 1 (a;) data by 
inversion. 



nounced with increasing temperature an has disappeared 
at T = 200 K. At T = 295 K the I 2 x{u) spectral density 
takes on the form of a simple MMP-form as was proposed 
by Millis el al. 36 for a spin-fluctuation spectrum with 
its maximum at ~ 80 meV. In these respects the 7 2 x(w) 
spectral density found for the optimally doped PCCO 
single crystal follows the pattern found for all other hole 
doped High-T c superconductors analyzed so far. It is, fi- 
nally, important to point out that the T = 5 K shows an 
additional peak centered around ~ 12.5 meV. This peak 
describes a possible coupling of the charge carriers to a 
mode at this energy. 

We add Figs. 5 and 6 to make two important points. 
Fig. 5 shows the t~^(uj) data for T = 30 K (solid line) for 
the single crystal. Superimposed are two straight lines. 
The first, (red) dashed line starts at the u> = offset 
at the value of the residual scattering rate t~^(uj = 0) 
and follows experiment up to about 45 meV. The second, 
(red) solid line starts near the origin and follows the data 
for energies > 45 meV emphasizing the change in slope 
at uj — 45 meV. If the first were perfectly flat (which it is 
not) the two straight lines would represent two processes, 
the Drude (coherent) part of the scattering rate and the 
boson assisted (incoherent) part, respectively. If for the 
latter we assumed that the boson is an Einstein mode 
at some energy uje which is large enough that the bo- 
son assisted processes are well separated from the coher- 
ent Drude part, then the coherent part of t~^(uj) would 
be constant and equal to the residual scattering rate at 
uj = until uj = uje is reached. At uje the boson assisted 
absorption sets in as an additional process and t~^(uj) 
develops a kink as is suggested by the intercept of the 
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FIG. 6: (Color online) The optical scattering rate t~^(lu) vs 
uj for the optimally doped (x — 0.15) PCCO single crystal 
in the superconducting state at T = 5K. The (black) heavy 
solid line corresponds to experiment. The two (black) light 
solid lines are used to emphasize the change in slope in r^p (ui) 
at uj ~ 17meV. The (red) dash-double dotted line represents 
the theoretical Top (uj) calculated using the spectral density 
I 2 x(v,T — 5K) (blue, dash-dotted line) for unitary impu- 
rity scattering with T + = 3.1 meV and c = 0. The (green) 
dotted line shows the theoretical r ( ^, 1 (w) calculated using the 
spectral density I 2 \{v, mod) (blue, dashed line) for the same 
unitary impurity scattering. The v-satis was displaced by the 
amplitude A = 5 meV of the superconducting gap at T — 5 K 
with respect to the w-axis. 



dashed line). As required by theory, the z/-scale was dis- 
placed by the amplitude of the superconducting gap at 
T = 5K, A = 5meV. Following the argument already 
applied in our discussion of Fig. 5 it becomes obvious 
that the peak at v = 12.5 meV in the spectral density 
I 2 x{v,T = 5K) coincides with the change of slope of 
the optical scattering rate at ~ 17meV as indicated by 
the intersection of the two (black) thin solid lines. Thus, 
this peak represents, indeed, a real feature of the optical 
scattering rate and describes the coupling of the charge 
carriers to a (weak) mode at 12.5 meV. 

We added two more curves to this figure. The (red) 
dash-double dotted curve represents the solution of the 
full non- linear Eliashberg equations (A.l) found for the 
spectral density I 2 x(v,T = 5K) in the superconducting 
state at T = 5K. The impurity parameters were set to 
r + = 3.1 meV and c = 0. [It is identical to the (black) 
thin solid line in Fig. 4(a).] We see that this result fol- 
lows closely the data over the whole energy region. The 
second, (green) dotted curve is also for the supercon- 
ducting state at T = 5K but now the spectral density 
I 2 x{v, mod) has been used. The impurity parameters 
remain unchanged. This results falls well below experi- 
ment in the region 17 < uj < 40 meV, thus demonstrating 
the importance of the 12.5 meV peak in I 2 x(v, T = 5 K) 
for best possible data reconstruction. It is worth not- 
ing that the dashed-double dotted and the dashed curves 
meet only at higher energies around 70 meV. This demon- 
strates clearly that a peak in I 2 x{uj) centered around the 
energy lue not only produces a kink in (uj = uje) but 
is also responsible for increased scattering at all energies 
above uje- 



two straight lines in Fig. 5. This corresponds precisely to 
the peak at the same energy seen in the I 2 x(uj) spectrum 
shown as a (blue) dash-dotted line in Fig. 5 with the hor- 
izontal (black) light dotted line indicating the peak cen- 
ter. The kink in T~p(uj) so identified is a clear signature 
of the peak at 45meV in I 2 x(oj)- Obviously, the (red) 
dashed straight line has in our case a non-zero slope as a 
function of energy. This is the signature of additional bo- 
son assisted processes even at energies u < 45 meV and, 
consequently, the spectral function I 2 x{oj) has non-zero 
weight at low energies. Next, we will examine this region 
more closely. 

In Fig. 6 we address the question whether or not the 
low frequency peak at 12.5 meV in the low temperature 
I 2 x{v,T = 5K) spectral density (blue, dash-dotted line) 
is really essential for data reconstruction in the region 
< uj < 40meV and reflects features contained in the 
T~p(uj,T = 5K) data. The (black) heavy solid line rep- 
resents the data and, as in Fig. 5 we use two (black) 
thin straight lines to emphasize the change in slope which 
takes place at ~ 17 meV in r^ 1 (w,T = 5K). For compar- 
ison we also present the spectral density I 2 x{ v , mod) in 
which the peak around 12.5 meV in spectrum I 2 x{ v i T — 
5 K) has been replaced by a straight line connecting the 
origin with the bottom of the valley at 25meV (blue, 



B. The x = 0.15 thin PCCO film 

We followed the procedure outlined in the previous 
subsection to invert the optical scattering rate t~^(uj) of 
an optimally doped, epitaxially grown thin PCCO film 
(x = 0.15, T c = 21 K) described in Sec. II B. The best 
over all data reproduction of the normal state data re- 
quired an impurity parameter t + = 5.0 meV which cor- 
responds to an impurity scattering rate t^ p = 31.4 meV 
(~ 253 cm -1 ) in excellent agreement with the experimen- 
tal value of 32.5 meV (~ 262 cm- 1 ) at T — 25 K. In the 
superconducting state the best agreement was found for 
an intermediate impurity scattering strength described 
by the parameters T + — 7.8125 meV and c = 0.75. We 
also note that in this case the clean limit critical tem- 
perature T c q « 57 K. This is a remarkable result as it 
compares well with the intrinsic critical temperature of 
61 K found for the single crystal. Thus, we conclude that, 
indeed, both samples present the same physical proper- 
ties typical for optimally doped PCCO. 

In discussing the results of our data analysis we pri- 
marily concentrate on the similarities in the electron- 
boson spectral densities I 2 x(uj) for the two optimally 
doped samples investigated here. For this purpose Fig. 7 
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FIG. 7: The electron spectral density I 2 x(^) vs u> for the 
PCCO x = 0.15 single crystal (solid lines) and thin film 
(dashed lines). The 7 2 x(^) of the thin film has been scaled 
by a factor of 0.5 while the single crystal data is unchanged. 
This emphasizes similarities between the two sets of spectra. 
The results are (a) for the superconducting state at T = 5 K, 

(b) for the normal state at T = 30 K in the single crystal is 
compared with the T = 25 K spectrum for the thin film, and 

(c) for the normal state at T = 100 K. 



presents the J 2 x(w) spectra for the x = 0.15 thin film 
(dashed lines) scaled by a constant factor of 0.5 while the 
corresponding spectra of the single crystal (solid lines) 
stay unchanged. Figure 7(a) is for T = 5 K in the super- 
conducting state and the two spectra reveal very similar 
features as functions of to. There is a first resonance 
like peak at <~ 14meV (~ 12meV in the single crys- 
tal) which corresponds in energy to the spin one reso- 
nance reported by Wilson et al. 23 and to the peak at the 
same energy found by STM experiments 25 in the sys- 
tem PLCCO. There is a second, main peak at 50meV in 
the thin film sample. The single crystal J 2 x(w) spectrum 
shows this peak at the slightly lower energy of ~ 45 meV. 
There exists, as yet, no experimental evidence from neu- 
tron scattering for the existence of a corresponding peak 
in the imaginary part of the spin susceptibility at such an 



energy in PCCO and, therefore, our analysis represents 
a specific prediction. This second peak is followed by a 
valley-hump feature which starts at ~ 80mcV with the 
hump centered at ~ 107 meV. After this structure the 
spectral density levels off to a weak background which 
extends beyond to = 300 meV and which has little struc- 
ture in the case of the thin film. This result is very similar 
to what has been reported recently for the system LSCO 
by analysis of optical data 21 (including the valley- hump 
structure) and by neutron scattering. 19 It is interesting 
to note an additional valley-peak structure beginning at 
~ 55 meV in the single crystal spectrum which cannot be 
observed at higher temperatures. This structure is most 
likely due to not completely removed structures in the 
data around this energy. 

Figure 7(b) shows the single crystal results for T = 
30 K (solid line) and the thin film result for T = 25 K 
(dashed line) in the normal state; both scaled as de- 
scribed above. While at T = 25 K a small signature 
of the 14 meV peak remains, it has basically vanished at 
T = 30 K in the single crystal (only a small shoulder 
is left). These differences could also be due to the very 
different nature of the two samples. The second peak 
is still very pronounced and stays basically unchanged 
in energy in both samples. The valley-hump structure 
which follows this second peak at around 80 meV stays 
unchanged in comparison to the T = 5 K spectra as is 
the background for lu > 200 meV. 

Finally, Fig. 7(c) presents the rescaled spectra for 
T = 100 K. We see that all structures have been smeared 
out in the spectra and the positions of the main peak 
change only little. The valley-hump structure which fol- 
lows the main peak has been smeared out almost com- 
pletely in the J 2 x(w) spectral density of the single crystal 
while it is still well developed in the thin film spectral 
density. It is also remarkable that the thin film spectra 
have a largely suppressed background for lu > 180 meV in 
comparison to the single crystal spectra. Increasing the 
temperature further results in spectra which are MMP- 
form 36 like and the maximum at T = 300 K moves to 
<~ 80 meV in the thin film and to ~ 90 meV in the single 
crystal. This corresponds to what has been observed in 
optimally doped Bi2212. 18 

We also note that the maximum of the hump in the 
T = 5 K and 25 K thin film spectra is at lower energies 
as compared to the same structure in the single crys- 
tal spectra. At 100 K, though, they agree, i.e.: in the 
x = 0.15 thin film spectra the maximum in the high 
energy valley-hump structure moves towards higher en- 
ergies with increasing temperatures while it stays fixed 
in the single crystal. 



IV. SUMMARY 

Motivated by the report of a spin one resonance in 
the imaginary part of the spin susceptibility in the 
electron-doped system PLCCO at the very low energy of 
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~ 11 meV by neutron scattering we studied the optical 
scattering rate reported for optimally doped PCCO sam- 
ples using the maximum entropy technique to extract in- 
formation on the electron-boson spectral density 7 2 x(w). 
This spectral density contains information on the cou- 
pling of the charge carriers to bosonic modes. We found 
that in the superconducting state the electrons couple 
to a bosonic mode centered around 12(14) meV and that 
there is a second higher energy group of modes centered 
around 45(50) meV which has not yet been observed by 
neutron scattering. Above T c the thin film sample still 
shows some coupling to the low energy mode at 25 K 
and while in the single crystal no coupling to this mode 
can be observed at T — 30 K which can certainly be 
due to the differences in the two samples. The second, 
high energy peak is clearly developed above T c in all 
samples and evolves into an MMP-form like spin fluc- 
tuation background with further increasing temperature. 
All I 2 x(co) spectra extend to energies > 300 meV. These 
results resemble very closely to what has been observed 
in optimally doped LSCO samples and, apart from the 
low energy mode, what has been reported for optimally 
Bi2212 samples. All this proves that the electron-doped 
system PCCO behaves in its charge dynamics much like 
all the other hole-doped High-T c cuprates. There is one 
important difference though, PCCO in contrast to all 
hole-doped High-T c superconductors investigated so far, 
is in the dirty limit. The residual scattering in all sam- 
ples investigated here is sufficiently large to substantially 
reduce the value of the critical temperature over its pure 



(intrinsic) limit. This fact is in full agreement with re- 
sults reported by Dagan et al. 37 from their analysis of 
PCCO-lead tunneling junctions. 
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APPENDIX: d-WAVE ELIASHBERG EQUATIONS 
FOR IMPURE SYSTEMS 

The generalization to a d-wave gap has already been 
published by Jiang et al. 38 and has been used to de- 
scribe various aspects of the superconducting state in 
the cuprates. In the mixed representation of Marsiglio 
et al. 39 they are of the form 



A(v + iO + ;$) = nTg cos(2tf) [X(u - iu m ) + \{u + iu m )} h+(iu m ) (A. la) 

m=0 

oo 

+ing J dz cos(2$)/ 2 x(z) [n(z) + f(z — v)\ h-{iw m — > v — z + i0 + ), 



and, in the renormalization channel, 

oo 

lo(v + i0 + ) = v + iirT ^2 [H p - iuj m) - A(zv + iu) m )] g + {iu) m ) 



m=0 



oo 

+iir J dzl 2 x(z) [n(z) + f(z - v)\ g-(icOm — ► v - z + M + ) 

— oo 

c 2 +g 2 _(ioj n -> v + iQ+) + h 2 _{iu n -> v + y ' 



Here 



/ A(iw m ;-!?)cos(2t?) \ / w(iw m ) 



I ^ / ' »=nv» w my \ I ^ 

\Ju) 2 (iuj m ) ± A 2 (iuj m ;i3)/ ^ \ y cu 2 (iuj m ) ± A 2 (iw m ;z9), 



r 



and the parameter g in Eq. (A. la) allows for a possi- channels. It is fixed to get the measured value of the 
ble difference in spectral density between the u> and A 
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critical temperature. In the above A(zu; m ;i?) is the pair- 
ing energy which its evaluated at the fermionic Matsub- 
ara frequencies u) m = nT(2m — l),m = 0, ±1, ±2, . . . 
and u>(ioj m ) are the renormalized frequencies evaluated 
at the same Matsubara frequencies; f(z) and n(z) are 
the Fermi and Bose distribution, respectively. Further- 
more, the d dependence of the pairing energy is de- 
scribed by A(iuj m ,$) — A(iu) m ) cos(2$) with d the po- 
lar angle in the two-dimensional CuO Brillouin zone. 
The brackets (•••)# are the angular average over and 
\{v) = J^dnaPFW/tu -fi + *0+). Eqs. (A.l) area 
set of nonlinear coupled equations for the renormalized 
pairing potential A(v + «0 + ;t?) and the normalized fre- 
quencies uj{v + i0 + ) with the gap A{v + i0 + ; r d) = A{v + 
i0 + ; r 9)jZ(y), where the renormalization function Z{y) 
was introduced in the usual way as uj(v + i0 + ) = vZ(y). 
Finally, A(iw„,i?) and u)(iui n ) are the solutions of the 
equivalent equations formulated on the imaginary axis. 40 
Impurity scattering is described by the term propor- 
tional to T+ in Eq. (A. lb) and enters only this equation 



because we assume a pure <i-wave model for the pairing 
potential with zero average over the Fermi surface while 
the impurity scattering is assumed to be isotropic. Here, 
T+ is proportional to the impurity concentration and c 
is related to the electron phase shift for scattering off the 
impurity. For unitary or resonant scattering c is equal 
to zero while c — ► oo gives the Born approximation, i.e.: 
the weak scattering limit. In this limit the entire im- 
purity term reduces to the form iirt + g(iu>„ — > v + i0 + ) 
with c absorbed into the impurity parameter t + . In 
the normal state A(u + i0 + ; r d) = and there is no 
need to distinguish any longer between unitary and Born 
limit impurity scattering. The scattering term reduces 
to iirt + sgn[Cj(iuj n —> v + i0 + )]- At the critical tem- 
perature linearized Eliashberg equations are valid, i.e.: 
h(iu) n ) ~ and g(iuj n ) — sgn[d)(iw n )]. Thus, at T c , 
t + {c 2 + sgn[a)(icj rl )]} = T + and this relates immediately 
r + to the impurity scattering rate via = 2irt + . 

The optical conductivity follows from knowledge of uj 
and A. The formula to be evaluated is 



2 i oo i 

a op {T,v) = ^UjdujtB^{^^[J{uj,v)-J{-uj,v)]\ . (A.2) 

The function J(co, v) is given by 

1 - N(u>; $)N(uj + v;d)- P(cu; d)P(cj + 1 + N*(uj; d)N(uj + y- g) + P*(cu; ti)P(cu + v; g) 



with E(w;0) = Ju> 2 (uj + i0+) - A 2 (u + iQ+;d), A(uj + iO+;-d)/E(uj;d). Finally, the star refers to the 
N{ U ;4) = *{» + iO + )/E(^), and P(u>;0) = c ° m P^ con ^te. 
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